Combined analysis of the integrated Sachs— Wolfe effect and cosmological implications 
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We present a global measurement of the integrated Sachs-Wolfe (ISW) effect obtained by cross- 
correlating all relevant large scale galaxy data sets with the cosmic microwave background radiation 
map provided by the Wilkinson Microwave Anisotropy Probe. With these measurements, the overall 
ISW signal is detected at the ~ 4.5(t level. We also examine the cosmological implications of these 
measurements, particularly the dark energy equation of state w, its sound speed c^, and the overall 
curvature of the Universe. The flat ACDM model is a good flt to the data and, assuming this model, 
we find that the ISW data constrain (1^ = 0.20j;Q jf at the 95% confidence level. When we combine 
our ISW results with the latest baryon oscillation and supernovae measurements, we find that the 
result is still consistent with a flat ACDM model with w = —1 out to redshifts z > 1. 



PACS numbers: 98.80.Es, 98.54. Aj, 98.70. Vc 

I. INTRODUCTION 

There is growing evidence that the expansion of the 
Universe is dominated by an unknown dark energy, which 
accounts for more than 70% of the total matter that we 
observe. Measurements of the Hubble diagram of type 
la supernovae [l|, Q indicate this dark energy is caus- 
ing the Universe's expansion to accelerate, while cosmic 
microwave background (CMB) anisotropics power spec- 
trum experiments, such as WMAP 3], limit the amount 
of spatial curvature which might otherwise explain the 
low dark matter density. Understanding the nature of 
this dark energy is essential, and many possible explana- 
tions have been proposed, from a cosmological constant, 
to quintessence, to modifying the laws of gravity. 

To discriminate between possible explanations of dark 
energy, many alternative probes have been developed, in- 
cluding measuring the Universe's geometry through ob- 
servations of the baryon oscillation (BAO) scale [j, [^, @| 
and probing the dark matter density through observa- 
tions of clusters of galaxies 0, @]. Here we consider 
another impact of dark energy, the creation of CMB 
anisotropics at late times via the integrated Sachs- Wolfe 
(ISW) effect [ol. While most CMB anisotropics are gen- 
erated at very early times, further fluctuations can be 
induced gravitationally at late times as photons pass 
through evolving gravitational potentials. If dark mat- 
ter dominates, the gravitational potentials do not vary 
with time, but the presence of dark energy or spatial cur- 
vature will cause the potentials to evolve at late times, 
producing new temperature fluctuations at low redshifts 
(primarily at z < 2). 

Directly observing these new CMB temperature 
anisotropics is challenging, primarily because their am- 
plitudes are a fraction of the anisotropics arising from 
higher redshifts. They are also predominantly seen on 



large scales, where the uncertainty in the observations 
from cosmic variance is biggest. The search for the ISW 
effect has instead focused on finding correlations between 
the CMB temperature maps and maps of the density 
which trace the local gravitational potentials jlOj. Pri- 
mordial anisotropies should be uncorrelated with the lo- 
cal density, making it possible to pull out the weaker ISW 
anisotropies. 

Many groups have detected this correlation using the 
accurate WMAP CMB data and various density probes 
distributed at a range of redshifts and in different regions 
of the electromagnetic spectrum: from shallow infrared 
observations to optical surveys such as the APM 

and the SPSS jl3l Jl4lT5l.ll^[T7t . radio galaxy catalogues 
[H, [3 , X-ray survevsJlSj. and the deepest quasars from 
the SDSS [23|. These measures span a range of redshift 
going from z = 0.1 to z = 1.5, where the ISW effect 
has been measured at significance levels typically around 
2 — 3cr, and appear generally compatible with the expec- 
tation from the ACDM model. 

Although indicative of the presence of dark energy, 
none of these measures alone has significant power to 
constrain models due to their low significance. Thus, it 
is important to combine the various observations; but 
some care must be taken in doing so. The surveys are 
often overlapping both in sky coverage and in redshift 
range, meaning there are likely covariances between them 
that may be important when considering a large scale ef- 
fect like the ISW. In addition, these measurements have 
been made with a variety of techniques, using angular 
correlations, Fourier modes, or a range of wavelet tech- 
niques [U, 122, 1131 ■ The error bars themselves have 
also been estimated using different techniques, using both 
jack-knife approaches and Monte Carlo simulations of the 
CMB sky. 

A combined analysis has been attempted in the past 
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adding several measures in order to extend the con- 
straining power in redshift and learn more about the 
behaviour of dark energy and other cosmological param- 
eters [13, HI]. However, this analysis largely ig- 
nored the differences in the observations and accounted 
for the covariances between experiments in a fairly ar- 
bitrary way. Here we perform a combined analysis by 
reanalysing all the observations in a consistent way, mea- 
suring directly the covariances between the different ob- 
servations using a number of different methods and look- 
ing at the cross-correlations between all the various data 
sets. In this way we hope to give a definitive result for 
the ISW evidence for dark energy and the resulting cos- 
mological constraints. 

This paper is structured as follows: we begin in sec- 
tion ^Jl] by giving a brief theoretical description of the 
ISW effect, how the cross-correlation is measured, and 
the important issue of estimating the covariance between 
observations. In section miU we describe the catalogues 
used for the cross-correlation, and in section i jlVl we show 
the measurements of the various cross-correlation func- 
tions between the different catalogues and their cross- 
correlation with the CMB. We discuss the significance 
of the measurements in section ^|V] and show the result- 
ing cosmological constraints in section Wli before some 
concluding remarks in section Will 



II. METHOD 
A. The ISW effect 

We begin with some theoretical background on the 
ISW effect and its detection. Most of the CMB 
anisotropics we observe were created at the redshift of 
last scattering {z ^ 1100) when the Universe was 400 
ky old, as a result of fluctuations in the photon den- 
sity, velocity and potential energy. Since that time, 
the CMB photons have travelled largely untouched, but 
anisotropics can still be produced gravitationally if the 
photons pass through time varying potential wells, 

^(n) = -2y'$[r,n(ro-T)]dr, (1) 

where r is the conformal time, the dot represents a con- 
formal time derivative and the integral is intended along 
the line of sight of the photon. (Throughout, we work in 
units where the speed of light is unity and for simplicity 
assume that the effects of anisotropic stress can be ig- 
nored). When the gravitational potential decays due to 
its linear evolution, this is usually referred to as the inte- 
grated Sachs- Wolfe effect, and if the potential decay is a 
result of non-linear evolution, as in clusters, it is referred 
to as the Rees-Sciama effect [29|]. 

The physical picture is very straightforward; as a CMB 
photon falls into a gravitational potential well, it gains 
energy; as the photon climbs out of a potential well, it 



loses energy. These effects exactly cancel if the poten- 
tial is time independent, but can result in a net kick if 
the potential evolves as the photon passes through it. In 
particular, we know that during the matter dominated 
era the gravitational potential stays constant, and so 
$ = 0, which means that in that era there will not be any 
ISW produced. This changes if dark energy or curvature 
become important at later times: in this case, <& 7^ 
and additional CMB anisotropics will be produced. In 
the usual case, the potential amplitudes decrease at late 
times, so that a temperature increase results from pass- 
ing through potential wells, while a temperature deficits 
results from traversing potential hills. 

Since we know from the CMB experiments such as 
WMAP that the Universe is very close to flat, we can 
attribute most of the late ISW to dark energy; observ- 
ing it can therefore provide a probe of dark energy, its 
properties and its evolution in time. In particular, the 
ISW effect could be one of the few ways of measuring 
the sound speed of dark energy [s^]. However, it can 
also provide useful information about the curvature of 
the Universe [3]| . 

Unfortunately, measurement of the ISW effect is made 
difficult by the embedding of the small ISW signal in 
the much larger (10 times) primary CMB anisotropics. 
Furthermore, the total ISW signal is due to all the den- 
sity fluctuations, both positive and negative, along the 
line of sight. On small scales, the individual tempera- 
ture differences are small and they tend to cancel out. 
The most significant ISW effect results from the coher- 
ent large scale potentials, but unfortunately these scales 
are precisely where cosmic variance is most troublesome. 

This problem can be overcome by examining how the 
ISW temperature correlates with the density of galaxies, 
which should trace the potential wells and hills which 
bring about the anisotropics. The observed galaxy den- 
sity contrast in a given direction ni will be 

Sg{hi) = J bg{z)^{z)6„,{ni,z)dz, (2) 

where dN / dz is the selection function of the survey, hg{z) 
its galaxy bias relating the visible matter distribution to 
the underlying dark matter, and 5m the matter density 
perturbations. Since the density (5,„ is related to the po- 
tential $ by the Poisson equation, the observed galaxy 
density will be correlated with the ISW temperature fluc- 
tuation in the nearby direction 112, which is 

^(n2) = -2 / e-<^)^{h^,z)dz, (3) 

where e~'^^^^ is the visibility function of the photons, 
which accounts for the effect of photons re-scattering fol- 
lowing reionisation. 

The galaxy bias, bg{z), can evolve in time or as a func- 
tion of scale; however, we will generally assume that it is 
time and scale independent for simplicity. For our pur- 
poses, a time dependent bias is equivalent to changing 
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the selection function of the survey. Scale dependence of 
the bias is more problematic, but on the very large scales 
(> 10 Mpc) we are considering, the scale dependence is 
expected to be weak 32, 33]. 

Given a map of the CMB and a survey of galaxies, the 
angular auto-correlation and cross-correlation functions 
are defined as 
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(4) 
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with the average carried over all the pairs at the same 
angular distance = |ni — n2|. 

It is possible to express these quantities in the har- 
monic space with the use of the Legendre polynomials 

Pr- 
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and the auto- and cross-correlation power spectra are 
given by 

= 4./ f A^(fc)/P^(fc)/f(fc) (7) 

(8) 
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4^ / -A2(fc)/f(fc)/f(fc), 



where A(fc) is the scale invariant matter power spectrum 
A2(fc) = 47rfc3p(fc)/(27r)3 and the two integrands are re- 
spectively 
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bg{z) — {z)6^{k,z)ji[kx{z)]dz, (10) 



where ^k, Smik, z) are the Fourier components of the 
gravitational potential and matter perturbations, ji{x) 
are the spherical Bessel functions and x is the comoving 
distance. These are calculated using a modified version 
of cmbfast SO. 



B. Theoretical signal-to- noise 

Unfortunately, the ability to detect the cross- 
correlation is limited because the signal falls off on small 
scales. Not only is cosmic variance an important factor, 
but there is also the problem of accidental correlations 
between the galaxy surveys and the CMB anisotropics 
produced at last scattering. Many independent measure- 
ments are needed to reduce the impact of such accidental 
correlations. The signal to noise ratio of the CCF with a 
particular survey is given by 
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For the ISW, we are usually in the weak correlation 
regime, so that Cj^UlcfCj'^ < 1. 



The signal to noise can be separated to obtain the 
contribution as a function of redshift; for a standard 
ACDM cosmology, most of the signal is expected at 
2; < 3, with the peak around a redshift of z ~ 0.5 [s^ . 
While the signal is highest at low redshifts, more inde- 
pendent volumes are available for higher z. The signal 
to noise scales roughly as the square root of the fraction 
of the sky observed. 

The most optimistic case is when the distribution of 
galaxies follows precisely the evolution of the ISW ef- 
fect. In this case, = C^^ = Cl^"^ where C/^^ is 
the spectrum of the ISW temperature anisotropics alone, 
which is assumed to be much smaller than the total CMB 
anisotropy, G^^^ <^ Cf^. Thus, the signal to noise re- 
duces to [lOl 
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This gives an optimistic total SjN ~ 7—10 for a standard 
ACDM cosmology. The ISW constraints which might 
arise from realistic future surveys can be found in [3a |. 



C. Correlation estimators 

Our aim is to estimate the correlation between sev- 
eral galaxy surveys and the CMB: as described above, 
this measure can be performed in the real space us- 
ing the CCF or in the harmonic space with the cross- 
correlation power spectrum. The two methods are theo- 
retically equivalent for a full sky analysis and both have 
been used to detect the ISW cross-correlations. However, 
when one moves away from the ideal full-sky scenario, it 
is more straightforward to account for the sky mask using 
the real space correlations, and therefore we will follow 
this approach here. 

The matter density and CMB temperature as well as 
their projections onto the celestial sphere are in principle 
continuous fields; however, we only have access to the 
sampling of these fields experimentally obtained by mea- 
suring the CMB temperature in some fixed directions and 
counting the number of galaxies in a given patch of sky. 
In practise, we pixelise these maps using the HEALPix 
pixelisation scheme [s^] , using a relatively coarse resolu- 
tion: iVsidc — 64, corresponding to A^pix = 49,152 pixels 
with dimensions 0.92° x 0.92°. This resolution is sufficient 
for a large scale correlations like the ISW and makes it 
tractable to perform large numbers of Monte Carlo sim- 
ulations. A finer resolution (iVgido = 128) was explored, 
but the results did not change significantly. 

In making the maps, we assign the average tempera- 
ture or the total number of galaxies to each pixel. The 
maps are masked according to the particular require- 
ments for each catalogue and the most relevant fore- 
grounds as discussed below. It is inevitable that some 
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pixels are only partially covered in the original survey, 
either because only part of the area was observed, or be- 
cause some of this area was masked out. In such cases, 
predominately occurring on the edge of the survey, the 
number of galaxies in a pixel is estimated as = n"^^ / fi 
where fi is the fraction of the pixel observed. The mask 
was obtained through sampling all objects in each cata- 
logue in a higher resolution (iVhigh = 512) as described 
in [20|. 

From these maps, both the auto- and cross-correlations 
were estimated, down- weighting those pixels with partial 
coverage proportionally to fi. For the auto-correlation 
functions (ACFs), we used the estimator. 
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fi (ni - n) fj {uj - n) , 



(13) 
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where n is the average number of galaxies in a pixel for 
the survey of interest, and = fifj is the weighted 
number of pairs of pixels with separation ??. For the 
temperature maps, we simply replace and n with the 
pixel temperature and average temperature of the CMB 
maps. 

More generally, we are interested in the cross- 
correlation function between the survey p and the sur- 
vey q; this is estimated similarly, accounting for the fact 
that the pixel weighting and mean number per pixel will 
depend on the survey. 



1 



(14) 



The number of pairs of pixels at a given separation, 
-^1' — J2ijfifj^ ■^ill depend on both of the surveys 
under consideration. This again extends to the density- 
CMB CCFs in the obvious way. 

We use Nf, — 13 angular bins in the range 0° < d < 
12°. We use a linear binning, and have explored the de- 
pendence of our results on the choice of binning, changing 
both the number and trying a logarithmic binning; nei- 
ther had significant impact on the results. 



D. Covariance estimators 

An important aspect of this calculation is the estima- 
tion of the covariance of the cross-correlation measure- 
ments. As described most recently by [l^, there are a 
number of different ways to calculate the errors on this 
measurement, each with their own advantages and draw- 
backs. Here, we calculate our errors in three ways: a 
Monte Carlo method (MCI), where the covariance ma- 
trix is estimated by measuring the CCF between random 
CMB maps while keeping fixed the observed density map; 
a second Monte Carlo method (MC2), similar to MCI 
but including also random density maps which are corre- 
lated at the expected level with the random temperature 
maps; and jack-knife errors (JK) which are estimated by 



looking at the variance of the CCF when patches of the 
sky are removed. 

The first approach is to generate random Monte Carlo 
maps of the CMB assuming the WMAP best fit cos- 
mology, and estimating the covariance matrix cross- 
correlating these maps with the true density maps 
(MCI). 

The WMAP third year fiducial model we use through- 
out this paper has baryon density il,b = 0.04185, matter 
density = 0.2402, Hubble constant Hq = 73.0, scalar 
spectral index = 0.958, optical depth t = 0.092 and 
amplitude of density fluctuations A = 0.80 at fc = 0.002 
Mpc-^ 

The MCI is the most widely used estimator in the lit- 
erature, though here we extend the usual calculation to 
account for covariances between the CCFs of the CMB 
with different surveys. This method is reasonably fast 
to implement and accounts for the cosmic variance and 
the accidental correlations with the CMB which are the 
primary source of error. However, it is asymmetrical, in 
that it does not account for the variance in the density 
maps or its Poisson noise; the MCI method also assumes 
there are no cross-correlations, though the expected (and 
observed) weakness of the cross-correlation indicate that 
this should not introduce a large bias. Finally, like all 
Monte Carlo approaches, it is model dependent and could 
fail if the data model is poorly understood (e.g. fore- 
grounds or non-Gaussianity of the maps). 

However, some of these problems can be avoided if we 
also generate random density maps for each catalogue 
based on the WMAP cosmology and the redshift distri- 
bution, with the addition of Poisson noise to the maps 
(MC2). In this case, we have the ability to account for 
the expected correlations between the maps as described 
in Appendix]^ This method is more time demanding, in 
that it requires more random maps for each correlation 
measurement; it also retains the unwanted model depen- 
dence, and unlike the previous method has no explicit 
dependence on any of the observed maps. 

To estimate the covariance between the different an- 
gular bins of a single CCF following the MCI and MC2 
methods for each catalogue k we use the following esti- 
mator of the full covariance matrix: 
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(15) 

where C'^^{di) are the mean correlation functions in the 
j-th angular bin over M realisations; the diagonal part 
of these matrices gives the variance of the CCF in each 
bin, Cfj — af, while the off-diagonal part represents the 
covariance between the points. 

The last method to estimate the covariance (jack-knife) 
consists in estimating the variance by generating mock 
density maps from the true ones, simply discarding a 
small patch of them. In practise, we can divide the orig- 
inal density map in M patches which have roughly equal 
area, and discard in turn a different patch to calculate 
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the CCF. The estimator for the covariance matrix is in 
this case, 
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X [cl'{^j)^C^'{^,)]. (16) 



The advantage of this method is its model independence, 
but it has the big drawback of giving different answers 
depending on the size and number of the discarded areas. 
It also implicitly assumes independence of the various 
patches, which is not always the case. 

Our ultimate goal is to measure the total covariance 
matrix between all the catalogues. To do so, we need to 
estimate the total covariance matrix C^', as the matrix 
that has in the diagonal blocks the single catalogue C^J , 
and in the off-diagonal parts is 
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(17) 

For simplicity, we redefine the indexes i,j in a way that 
they run from 1 to iVtot — -/Vbin x -^catj i-e. redefining the 
data, theory and mock arrays as the concatenation of all 
catalogues' CCFs with the CMB. In this way, the covari- 
ance matrix is simply the square matrix Cy , identical to 
the Eq. (jl7p but now with dimension A^tot- A similar 
expression can be easily obtained for the JK case. 



III. THE CATALOGUES 

To best detect the ISW effect through the cross- 
correlation technique, we ideally require surveys covering 
large fractions of the sky, so that accidental correlations 
will cancel out. The surveys also need to be sufficiently 
deep, in order to probe the gravitational potentials where 
the ISW effect is being created. Ideally, we would like 
to span the redshift range < z < 3, separated into 
subsamples of different depths so as to measure the red- 
shift dependence of the effect and get some handle on 
the evolution of the dark energy. However, only rather 
coarse redshift information is required, so redshift errors 
of Az ^ 0.1 obtainable through photometric methods 
should be sufficient for these purposes. This is beyond 
the present state of the observations, but the differences 
in the redshift distributions of the various samples does 
provide some limited information on the dark energy evo- 
lution. 

At present, the best surveys available for this pur- 
pose (and where ISW detections have previously been 
found) include the following: the optical Sloan Digital 
Sky Survey (SDSS), the infrared 2 Micron All-Sky Survey 
(2MASS), the X-ray catalogue from the High Energy As- 
trophysical Observatory (HE AO) and radio galaxy cat- 
alogue from the NRAO VLA Sky Survey (NVSS). The 



high quality of the SDSS data allows us to extract some 
further subsamples from it, consisting of Luminous Red 
Galaxies (LRG) and quasars (QSO) in addition to the 
main galaxy sample |39|. These are the samples we use 
in our analysis below, and include most of the significant 
reports of the ISW detection. Because the data has not 
been publicly released and since it is not significantly 
deeper than 2MASS, we omit the APM galaxy survey, 
which has also been reported to have evidence for ISW 
cross-correlations T^]. 

We show in Fig. [T] the redshift distributions dN/ dz of 
the catalogues we use, normalised to unity; we can see 
that they span a redshift range < z < 2.5, similar 
to the theoretical requirement, although the overlap be- 
tween different samples is significant. This means that 
the covariance between the measures could be large: one 
of the goals of this paper is to quantify it. 
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FIG. 1: The redshift distributions of all catalogues dN/dz 
normalised to unity. The significant overlap between redshift 
distributions (especially for the X-ray and radio surveys) re- 
sults in a covariance matrix with significant non-diagonal el- 
ements. 



In the rest of this section, we will present the charac- 
teristics of all the samples we use, in order of increasing 
redshift. 



A. 2MASS 

The 2 Micron AU Sky Survey (2MASS) is an infrared 
catalogue; its extended source catalogue (XSC) [iO] con- 
tains ~ 800,000 galaxies with median redshift z ~ 0.1 
and, unlike the point source catalogue (PSC) is almost 
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free of stellar contamination. Some evidence for ISW 
cross-correlations has been seen in 2MASS previously 
[ill . and we largely follow the galaxy selection of 
those previous analyses here. 

Accordingly, we select galaxies according to their Kg- 
band isophotal magnitude K2Q (k_m_i_20_c, - 20 mag / 
arcsec^). These magnitudes are corrected for Galactic 
extinction using the infrared reddening maps by [4l|. 
as ~ K2Q — Ak, where the extinction is Ak = 

0.367(i? — V). The requirement of completeness of the 
catalogue is satisfied by imposing a cut in magnitude 
< 14.0, while we can exclude low redshift sources 
with the condition > 12.0. We only include objects 
with a uniform detection threshold (use_src = 1), and 
remove known artifacts (cc_f lag 7^ a and cc_f lag ^ z); 
we also exclude a small fraction of objects where the mag- 
nitude or its error were not recorded. 

In addition to the pixelisation geometry mask, we fol- 
low earlier analyses 11, 12] excluding areas of the sky 
with high reddening, discarding pixels with Ak > 0.05; 
this leaves 69% of the sky and 718,000 galaxies after ex- 
cluding artifacts. It is reported by [111, [l3| that the red- 
shift distribution of these galaxies is well approximated 
by the function: 



dN 
dz 
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exp 
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(18) 



where the parameters are zq — 0.072, /3 = 1.752 and 
m = 1.901. This distribution is shown together with the 
others in Fig. [TJ 

To check the consistency of the dataset and its bias we 
calculate its auto-correlation function (ACF). The mea- 
sure is in good agreement with the predictions for the 
best fit WMAP model with a galactic bias bg — 1.4 as 
found by JJi] , as we can see in Fig. O 



B. SDSS galaxies 

The SDSS Sixth Data Release (DR6) [H, \^ is 
the largest wide optical galaxy survey available at the 
present for the northern hemisphere. From this cata- 
logue we select a magnitude limited subsample 18 < 
r* < 21; this catalogue contains 30 million galaxies. 
Here r* is the extinction corrected r SDSS iibercalibrated 
model magnitude, i.e. using the SDSS variables r* — 
ubercal.modelMag_r — extinctionjr: this corresponds 
to the procedure of fis'l, with the difference of using the 
sixth data release and the iibercalibrated model magni- 
tude instead of the Petrosian magnitude, which is less 
reliable for faint objects. We apply the pixelisation ge- 
ometry mask and, in addition, we discard the pixels most 
affected by reddening, with Ar > 0.18. We also discard 
the southern stripes, since they are most affected by fore- 
grounds and edge effects. 

We select only objects with photometric redshifts be- 
tween 0.1 < z < 0.9 and with an error on the redshift 



az < 0.5^, leaving 23.5 million galaxies in the catalogue. 
We could use these photometric redshifts as the basis 
of the theoretical calculations; however, since the distri- 
bution of the photometric redshifts can be affected by 
singularities in the redshift determination procedure, we 
use instead a fit to their distribution with the smooth 
function of Eq. (fT8|) . The best fit parameters are in this 
case Zq = 0.113, (3 — 1.197 and m = 3.457, corresponding 
to a median redshift Zmcd = 0.32. (The results are actu- 
ally independent of whether the fit or the actual redshift 
distribution is used). The fit is shown together with the 
others in Fig. [T] 

The ACF is in agreement with the prediction for the 
WMAP best fit cosmology and a bias 5g = 1, as we can 
see in Fig. [2l 



C. SDSS LRG 

Luminous Red Galaxies (LRGs) from the SDSS have 
been used often to find evidence for the ISW effect, as 
they have a deeper redshift distribution than the ordi- 
nary galaxies, with a mean redshift of z ~ 0.5 [isl ITgI ITtI . 
In this analysis we use the MegaZ LRG sample f4J, |45| 
which contains 1.5 million objects from the SDSS DR6 
selected with a neural network [4^] • To ensure complete- 
ness we require that i < 20. To reduce stellar contami- 
nation we implement cuts on Ssg, which is a variable of 
the MegaZ neural network estimator, defined such that 
Ssn = 1 if the object is a galaxy, and Ssg = if it is a star 
[ij]. Following the conservative suggestion by [ill, we 
choose a cut Sgg > 0.2, which is reported to reduce stel- 
lar contamination below 2% while keeping 99.9 % of the 
galaxies. Stricter cuts have been tried with no significant 
changes to the CCF. 

The mask we apply to this catalogue is a combination 
of the pixel geometry mask and two foreground masks, 
to account for seeing (cutting pixels with median seeing 
in the red band greater than 1.4 arcsec) and reddening 
(cutting pixels with median extinction in the red band 
Ar > 0.18). The redshift distribution function in this 
case is found directly from the photometric redshifts that 
are given in the catalogue, and is shown in Fig. [TJ 

We show the auto-correlation function in Fig. [21 where 
we can see that this is in agreement with the theoretical 
prediction from the best fit WMAP cosmology and a bias 
bg = 1.8, which is compatible with the estimate bg — 
1.7 ± 0.2 shown by [i^, although some excess power at 
large scales is present, which might be explained as being 
produced by a residual stellar contamination. 



D. NVSS 

The NRAO VLA Sky Survey (NVSS) is a fiux limited 
radio survey at a frequency of 1.4 GHz, with a mini- 
mum flux of ~ 2.5 mjy. It is complete for declinations 
S > —40°, covering roughly 80% of the sky and contains 
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1.8 • 10^ sources. The mask to this catalogue is a com- 
bination of the most aggressive WMAP mask (kpO) plus 
a cut around point sources as described in [43|, which 
also describes corrections made for a systematic in the 
mean density as a function of declination. The cross- 
correlations between NVSS and WMAP have been ob- 
served by a number of groups, both in the correlation 
functioi i fl8l.[l9l and using an array of wavelet techniques 
[2l|,[2l[2fl2i7 

The redshift distribution is uncertain; we base ours on 
models of Dunlop and Peacock [i^ which seem to be still 
widely accepted, and are largely consistent with obser- 
vations of cross-correlations with other surveys (though 
see below for further discussion) . We calculate the auto- 
correlation function and present it in Fig. [21 there is 
good agreement with the theory from the WMAP best 
fit model and a galactic bias 6o = 1.5, compatible with 
the result bg — 1.5 ± 0.2 by [43|, although we see some 
excess power at small scales. 



E. HEAO 

The High Energy Astrophysical Observatory (HEAOl- 
A2) data set is a full sky flux map of hard X-rays counts 
in the 2—10 keV energy range J49|. We use the map 
and the mask determined by |50l . |5l| : the map is masked 
for the galactic plane, a round area around the galactic 
centre and patch areas around bright point sources. The 
redshift distribution is also uncertain and provide d by 
modelling the X-ray background, as described in [50ll52||. 

The modelling of the theoretical ACF for this catalogue 
is more complex than those considered above, in that we 
are looking at flux rather than number counts and the 
experimental beam is large compared to the pixel size. 
(The point spread function of the beam is well modelled 
by a Gaussian with a full width, half maximum size of 
i?FWHM = 3.04° dll). In addition, the number of pho- 
tons is small, so there is an additional contribution from 
the photon shot noise. Thus, the observed correlation 
is the sum of three terms: the intrinsic correlations, the 
Poisson correlations due to finite numbers of sources and 
shot noise due to the finite number of photons. The vari- 
ance of the X-ray map is dominated by photon shot noise 
(41%) and Poisson correlations (45%) while intrinsic cor- 
relations are relatively small (14%). However, the shot 
noise contributes only to the 0° ACF while the Poisson 
correlations fall off more quickly with angle than intrin- 
sic correlations and become sub-dominant for 9 > 4°. 
Consequently, the combination of shot noise and Poisson 
correlations are not the primary component of the total 
noise in the ISW signal. We can see in Fig. [2] that the 
total modelled ACF fits the observations on large angles, 
assuming the WMAP best fit model and a galactic bias 
bg = 1.06, as found by ll8[. 



F. SDSS QSO 

The quasar survey we use comes from the SDSS DR6 
through the NBC-KDE catalogue by [Ullil, that con- 
tains over a million quasars. This new DR6 edition of 
the catalogue does not include as many parameter cuts 
as did the previous DR4 version. To obtain the cleanest 
possible dataset, and for consistency with our previous 
measure of the cross-correlation [23| , we only used quasar 
candidates selected via the UVX-only criteria used in the 
previous version of this photometric quasar catalogue. In 
addition, we consider only objects with a good (positive) 
quality flag. Following our previous results [20| we im- 
pose a cut in reddening, discarding areas with Ag > 0.18. 
After these cuts, we are left with N ~ 500, 000 quasars. 

This catalogue comes with estimated photometric red- 
shifts, upon which we base the redshift distribution 
shown in Fig. [T] There is evidence of some excess power 
in the ACF on large angular separations that indicate 
faint stars are still present in the catalogue after these 
cuts, as seen before in [2QI]. The amount of stellar con- 
tamination is ~ 3%, as found by [541] . from comparison 
with the ACF of a random sample of stars taken from the 
SDSS, and does not contribute to the correlation with the 
CMB, as expected. We can see in Fig.[5]the ACF for this 
sample; this is in good agreement with theoretical expec- 
tations and determines the bias of bg = 2.3, as previously 
found in [Hi^. 

IV. RESULTS 

In this section we present the measurements of the all 
the correlation functions between the data sets we con- 
sider and their covariance. 



A. Density-density cross-correlations 

We begin by examining the cross-correlations between 
the different density maps. These measurements are 
shown in Fig. [2l with the auto-correlation measurements 
along the diagonal. This is the first measurement of the 
cross-correlations between most of these data sets. The 
error bars are estimated by Monte Carlo realisations of 
all the data sets (MC2, as described above). 

The measurements largely agree with their theoreti- 
cal predictions, which are based on the WMAP best fit 
model using the visibility functions in Fig. [1] and a lin- 
ear bias for each. The agreement is to be expected for 
the auto-correlations, which were the basis for the esti- 
mates of the linear bias. However, the cross-correlation 
measurements provide a useful consistency check for our 
model, and in particular for the visibility functions, since 
the cross-correlations are most sensitive to the degree 
that the measurements overlap in redshift. 

The largest discrepancy between the measurements 
and theory is in the NVSS-2MASS cross-correlation. 
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FIG. 2: Measures of the two-point correlation functions between all the combinations of catalogues, where the units in the 
X-axis are degrees. The auto-correlations are on the diagonal, and the solid (red) lines show the theory from WMAP best fit 
cosmology and the galactic bias from the literature. The largest discrepancy with theory, in the NVSS-2MASS CCF, can be 
addressed by a small change in the assumed NVSS redshift distribution (blue dashed line). 



where the theory is roughly twice as large as expected. 
This is perhaps not unexpected, since the NVSS visibility 
function is known to be uncertain, and the overlap with 
2MASS is in a narrow region of redshift. It does indicate 
that less of the NVSS correlations are arising from the 
2MASS redshift range than expected in the model. This 
could be because either the low redshift tail of the NVSS 
visibility function is overestimated relative to the high 



redshift region, or because the bias of the radio galaxies 
increases as we move to higher redshift. This can be ad- 
dressed by a small change in the visibility function, as 
demonstrated by the blue dashed line in the panel (in 
this case we arbitrarily imposed a low redshift exponen- 
tial damping in the visibility function, leaving the rest 
unchanged). Such a change does not significantly affect 
the expected CMB cross-correlations considered here. 
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B. Temperature-density cross-correlations 

We next turn our attention to the cross-correlation 
functions between each density map and the CMB maps 
from WMAP 3. We use the internal linear combination 
(ILC) maps from WMAP , which are the cleanest data, 
although we have checked that the results do not depend 
on the frequency (see below), and we also apply the kpO 
mask to them, cutting the galactic plane region. As we 
can see in Fig. [31 the measures are again largely in agree- 
ment with the theoretical predictions for the WMAP best 
fit model. 

We have also checked the results obtained with the new 
WMAP 5 data, and we have not found any difference in 
the correlations. This is expected, since WMAP maps 
are already cosmic variance limited at large scales. 

We now discuss the results obtained following the three 
methods of error estimation discussed in iflll above. 

1. Temperature- only Monte Carlo errors 

We generate 5000 Monte Carlo simulations of the CMB 
anisotropy map with the WMAP best fit parameters. We 
estimate the covariance matrix for each catalogue using 
Eq. (lisp , and the total covariance matrix follows from its 
generalisation. 

These are the errors shown in the top panel of Fig. [31 
as we can see, the errors are quite large, especially for 
the low redshift catalogues, and the significance is further 
decreased by the high correlation between the points. We 
have checked that these errors converge; the convergence 
is already good after ~ 700 Monte Carlos for each single 
catalogue, and after ^ 3000 Monte Carlos for the full 
covariance matrix. The covariance between the points is 
shown in Fig. [51 

2. Full Monte Carlo errors 

In this case, in addition to 5000 new mock CMB maps, 
we also generate 5000 mock density maps for each cata- 
logue, correlated as expected theoretically, based on the 
WMAP best cosmology and their redshift distributions. 
In addition, the Poisson noise is added due to the ex- 
pected number of objects per pixel. 

The result calculated in this way is shown in the bot- 
tom panel of Fig. [31 and the relative full covariance ma- 
trix in the bottom panel of Fig. [51 We can see that the 
errors estimated in this way are generally consistent with 
their MCI counterparts. 

The largest difference between the approaches is in the 
covariance between the cross-correlations measured with 
diff'erent data sets (Fig. [5]). Using the observed density 
maps yields both positive and negative covariance, while 
the covariance is only positive when all the maps are sim- 
ulated. In the first approach, the strongest correlations 
are between the SDSS subsamples and 2MASS. In the 



second approach, which is purely theoretical, the largest 
covariances are between 2MASS, NVSS and HEAD. The 
NVSS-HEAO covariance is expected to be large, since 
they are both essentially all sky maps and have similar 
redshift coverage. The large covariance between 2MASS- 
HEAO and 2MASS-NVSS is more surprising given the 
differences in the redshift distributions, but seem to be 
driven by the low redshift tail of the NVSS and HEAD 
distributions. As noted above, the cross-correlations are 
smaller than expected theoretically for 2MASS-NVSS 
(and to a lesser extent for 2MASS-NVSS). This indicates 
that the overlap of 2MASS with NVSS and HEAO is less 
than assumed, and that we have likely overestimated the 
covariance somewhat. However, the low significance of 
the 2MASS CCF means this has a small impact on the 
final result. 

The differences between the two methods appear large 
for the off-diagonal elements. The reasons for these dif- 
ferences are unclear, but they suggest that the observed 
density maps are somewhat atypical of those simulated. 
However, it is not surprising that any particular realisa- 
tions would appear atypical in some way. Despite these 
differences, these covariance matrices give comparable fi- 
nal significance, as is discussed below. 



3. Jack-lcmfe errors 

For completeness, we also present the errors estimated 
from a jack-knife method. However, there is more am- 
biguity in implementing this method, leading to uncer- 
tainty in the resulting estimates in the errors. 

One issue is what patch size to use. Ideally for the 
jack-knife approach one would like the cross-correlation 
observations to be uncorrelated between patches. In re- 
ality, some correlation is inevitable. We also need enough 
patches to estimate the full covariance matrix without it 
becoming singular, which drives the size of the patches 
down. Thus, some kind of compromise is required. 

Since we are interested in the CCF on scales of a few 
degrees, we choose a patch size of order 10 square de- 
grees. Because the surveys have different geometry and 
masks, the number of sectors M will be different in each 
one. The number of patches we can have in this way is 
generally low (~ 100), so we cannot estimate the total 
covariance matrix which, having a dimension TVtot = 78, 
requires at least a few hundred independent random mea- 
sures to be correctly estimated. 

Cross-correlation measurements also introduce other 
issues, since the CMB and density maps are often cover- 
ing different regions of the sky. 

In the end, we tried to be conservative and ensure the 
most independence between the subsamples by only in- 
cluding data which were in the CMB and the density 
maps, and masking out both maps in the jack-knife es- 
timates. The results we obtained are shown in Table [H 
where we compare the results obtained with jack-knife of 
the density map only and of both density and tempera- 
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FIG. 3: Monte Carlo error estimation. Measurements of the cross-correlation functions between all the catalogues and the 
WMAP CMB maps (black points), compared with the theory from WMAP best fit cosmology and the galactic bias from the 
literature (red solid lines). The best fit amplitudes and their 1 — a deviations are shown in blue (dashed). In the top panel, the 
errors are calculated with 5000 temperature-only Monte Carlos and, in the bottom panel, Monte Carlos for temperature and 
density including expected correlations. We see that the errors are comparable for individual observations. Because of known 
contamination from the Sunyaev-Zeldovich effect in the 2MASS data [Tl[, the four smallest angle bins were excluded from the 
fits. 
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FIG. 4: Jack-knife error estimation. The lines are the same as in Fig. |3l The errors are somewhat smaller than seen from the 
Monte Carlo estimates, possibly due to correlations between the jack-knife subsamples. 



ture maps using identical masks. 

The jack-knife ambiguities are even more problematic 
when calculating the full covariance between observations 
using different density maps, since the density maps often 
will not overlap on the sky. For this reason, and because 
such a large number of jack-knifes are required to esti- 
mate the total covariance matrix, we do not attempt to 
estimate it here. 

The error bars we estimate using the jack-knife method 
are of the same order of magnitude as those seen in the 
Monte Carlo approaches, but are somewhat smaller lead- 
ing to higher significance in the detection. This could be 
due to the lack of independence of the jack-knife patches, 
or because some aspect of the data is missing from the 
Monte Carlo approach. We will use the Monte Carlo 
estimates below, focusing primarily on the results from 
MC2. 



C. Foregrounds & systematics 

Since the ISW effect is gravitational in origin, it is 
frequency independent as are the resulting CMB-density 
cross-correlations. However, a frequency dependence 
may in principle be introduced by foregrounds and lo- 
cal contamination, such as the SZ effect. In Fig. [5] we 



compare the CCF obtained with the different frequency 
bands from WMAP (ILC, W, V and Q bands), and we see 
that the result is substantially independent of frequency, 
with the exception of the 2MASS catalogue. However, 
the 2MASS CCF detection is of low significance and our 
final answers are not greatly sensitive to its inclusion. 

Foreground contamination of the ISW signal is gener- 
ally produced at low redshifts. A good way to make sure 
that such effects are not dominating the measurement is 
to check for the sensitivity to the masking of these fore- 
grounds (e.g. [111). For samples derived from the SDSS 
(galaxies, LRGs and QSO), we test for foreground effects 
by cutting the 20% of pixels with the highest redden- 
ing (extinction), seeing, sky brightness, and number of 
unresolved point sources. The most relevant masks are 
the reddening and seeing masks which do not substan- 
tially change the results. For the other samples (2MASS, 
HEAO and NVSS), we do not explore the masking, but 
we refer to the foreground analyses presented in earlier 
papers dH, [H, 113, iHH . 

D. Comparison with previous measures 

We briefly compare our CCF measurements to others 
in the literature. 
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FIG. 5: The total covariance matrix obtained with 5000 
Monte Carlos, normalised. The top panel shows the 
temperature-only Monte Carlos, while the bottom panel is the 
result of the full Monte Carlos. While the diagonal (single ex- 
periment) covariances are similar, those between experiments 
(off-diagonal) are somewhat different. 



nificant evidence for cross-correlations. 



2. SDSS galaxies 

The main galaxy sample from the SDSS has a mea- 
sured CCF which is also in good agreement with the the- 
ory. In this case, we note that we do not find agreement 
with the previous result of [l3| , who reported a measured 
CCF of almost double the amplitude that we detect. 

After discussions with the authors [l^, we jointly 
found this discrepancy resulted from an additional clean- 
ing cut, where they discarded all galaxies with a large 
error on their Petrosian r magnitude, imposing the con- 
dition petroMagErrjr < 0.2. Imposing this same condi- 
tion, we found that we could reproduce their result. Fur- 
ther, masking those areas with high proportion of Pet- 
rosian error also gave similar results. 

However, the motivation for such a cut is unclear. It 
is known that the Petrosian magnitudes are not accu- 
rate for faint objects, for which the best estimator is 
the model magnitude [1^. While having objects with 
a well measured magnitude is desirable, we see no reason 
why cutting galaxies on the basis of a poor estimate of 
their magnitudes should double the correlation with the 
CMB. This could happen if it were produced by some 
foreground mechanism, such as seeing or reddening, but 
we checked that none of the possible foreground maskings 
raised the CCF in any way comparable to the aforemen- 
tioned cut. 

Therefore, lacking a valid reason to include this cut, 
and preferring to be conservative, we do not make the 
Petrosian error cut and our CCF is thus lower than seen 
by Cabre et al. . While it is worrying that a choice of 
masking has such a dramatic effect on the amplitude of 
the observed cross-correlation, it should be noted that the 
cross-correlation was largely independent of other mask- 
ing choices. 



SDSS MegaZ LRGs 



1. 2MASS 

From Fig. [3] it is clear that the CCF for the 2MASS 
survey is consistent with zero. Previous analyses of 
these data found some evidence for a positive correla- 
tion [ill , [l^ : however, these were performed in Fourier 
space and included modelling of the SZ effect, which man- 
ifests itself with anti-correlations at small angular scales. 
Indeed, it appears in Fig. 3 that the observed CCF turns 
over at small angles. If the smallest four angular bins 
are removed, the fit to the CCF is consistent with the 
ACDM theory; however, it is only significant at the ~ la 
level. In any case, 2MASS appears to have the least sig- 



The result for the LRG is the highest in comparison 
with the ACDM theory. It agrees with the result of Ts*! . 
A direct comparison with [13] and [l6j is more difficult 
because these analyses use multiple photometric redshift 
bins. Concentrating on [l3| (since it also does its analy- 
sis in physical space, rather than Fourier space), we find 
approximately the same detection significance as their 
single redshift bin measurements for similar data sets. 
An updated version of this paper (available on the astro- 
ph archive, but also unpublished) calculates a global 
value using all four of their LRG samples, and detects 
a CCF with significance somewhat higher than we mea- 
sure in this work. This is likely due in part to a somewhat 
larger redshift baseline for their measurement as well as 
the fact that they calculated their covariance matrix us- 
ing a method similar to our MCI case. As one can see 
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FIG. 6; Comparison of the CCF functions obtained with the different WMAP frequency bands. The black (solid) is using the 
internal linear combination map; the blue (long-dashed) uses the W band, the green (short-dashed) uses the V-band and the 
cyan (dotted) uses the Q-band. The thick red curves show the ACDM prediction. 



from Fig. [5l samples which cover very similar areas and 
have significant redshift overlap (as is the case with their 
LRG photometric redshift samples) can result in stronger 
anti-correlation between samples than one observes in co- 
variance matrices generated with the MC2 method. This, 
in turn, would lead to a moderate over-estimation of the 
detection significance. 

4- Other measurements 

Not surprisingly, since we use the same maps generated 
from HE AO and NVSS, our results are in agreement with 
previous measures by [H, [l^, and the amplitudes are 
consistent with the theoretical predictions. As discussed 
above, the new Monte Carlo approach give consistent 
answers for individual experiments as the temperature- 
only Monte Carlo approach used in earlier analyses. 

We found that the measured CCF for the quasars is 
consistent with the earlier measurement and the expecta- 
tion from theory, and it is independent from the cleaning 
level of the catalogue. 

In conclusion, all the measured CCF agree with 
the previous results and with the ISW theory for 
a ACDM model, although they are in some cases 



marginally higher than theory predicts. 



V. SIGNIFICANCE OF THE RESULT 

Having established the measures of the CCFs and the 
total covariance matrix, we discuss the significance of this 
result and its consequences. 



A. Single catalogue significance 

Assuming that the detected cross-correlations are due 
to the integrated Sachs-Wolfe effect, we can assign a 
significance value to the measure if the errors on the 
cross-correlation are taken to be Gaussian. For each 
catalogue, we can compare the measured CMB-density 
cross-correlation C{'di) with the theoretical expectation 
obtained from the WMAP best fit cosmological parame- 
ters with our modified version of the cmbf ast code (341] . 

We perform the likelihood analysis first described 
in [50*1 . The shape of the CCF for each catalogue is 
assumed to follow the ACDM predictions. The theory 
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template is 



(19) 



where g{'di) is theoretical prediction of the WMAP best 
fit model and A is the fit amplitude, which will depend 
on the visibility function of the catalogue in question. 
Maximising the likelihood 



-1/2 



X exp 



,(20) 



we can find the best value for each A, 



A 



and the variance 



2^1.1=1 '-'11 9^9j 



N 



(21) 



(22) 



We can also simply obtain the signal to noise ratio as 
S/N = AjoA- 

The results obtained in this way with errors calculated 
with the three methods are summarised in Table HI and 
the resulting amplitudes and their errors can be seen in 
Fig. [3] and Fig. ID Here we have allowed a separate ampli- 
tude A for each catalogue. Note that while the observed 
CCF is the same for the different methods, differences 
in the covariance matrices can result in different best fit 
amplitudes. 

It is possible to check that the Monte Carlo estima- 
tion has converged after N realisations by estimating the 
uncertainty on the errors. In detail, we use a jack-knife 
approach consisting in observing the effect of removing 
M = 10 different subsets of the N = 5000 realisations of 
the MC2 method. The estimator of the uncertainty on 
SIN is 



M ■ 



M 



^ M 

-T.[^S/N\ 

i=l 



S/N 



(23) 



where (S/N)^ are the signal to noise ratios obtained with 
each subset of — M Monte Carlos, and S/N is their 
average. We find in this way that the uncertainty on the 
S/N is less than 5%, indicating the level to which our 
Monte Carlos have converged. 



-^tot = -^bin X Ncati running over each of the bins of the 
of the observed (and theoretical) cross-correlation func- 
tions for each of the density catalogues. Using the full 
covariance matrix, we can follow again the same proce- 
dure, and find a single best fit amplitude. 

The results obtained in this way are shown at the bot- 
tom of Table |T1 The significance of the two different 
Monte Carlo methods, MCI and MC2, are 4.4ct and 4.5cr 
respectively. We also find that the uncertainty on the 
S/N for the joint amplitude is again less than 5%. 

The two MC methods produce similar detection signif- 
icances, but this could be a lucky coincidence, since the 
covariance matrices relating different surveys are much 
different. Both methods suggest some pairs of observa- 
tions should be strongly correlated, but which pairs are 
strongly correlated depends on the method. If the covari- 
ance between surveys were ignored, the total significance 
would be about 5.8a. Perhaps it is not surprising then 
that adding similar levels of covariance between experi- 
ments with comparable individual detection levels would 
have a similar effect on the total significance. 

As in the case of fits to individual correlation functions, 
strong covariances can have results which are counter- 
intuitive. For example, the fit for the total amplitude 
using the MCI approach is smaller than any single sur- 
vey would suggest. Also, adding the small angle 2MASS 
CCF, believed to be suppressed by SZ, actually increases 
the fits by about 0.2a despite the points themselves be- 
ing lower than the theory. These effects suggest that 
the degree of covariance between the different measure- 
ments might be over-estimated, which would not be sur- 
prising given the much different systematics in each ex- 
periment. Even adding a small degree (5%) of diagonal 
noise is enough to increase the total MCI amplitude to 
A = 1.14±0.26, with a corresponding S/N = 4.4, so that 
it is more consistent with the amplitudes of the individual 
experiments. The MC2 result is not affected by such a 
change, because the total amplitude is already consistent 
with the individual survey measurements. 

Note that the theoretical model associated with a par- 
ticular best fit amplitude is not unique. While increasing 
the dark energy density will generally increase the ISW 
effect, the effect will generally be redshift dependent and 
could impact different catalogues differently. However, 
the ACDM model without any tweaking {A = 1) im- 
proves the likelihood at ~ 4.5tT compared to the absence 
of cross-correlations. Below we compare to specific alter- 
native cosmologies without any scaling amplitude. 



C. Tests 



B. Joint significance 

We can easily generalise this to combine the different 
catalogues and obtain a single significance. Redefining 
the indexes i,j in a way that they now run from 1 to 



Another way to assess the significance of the measure 
with respect to a theory is simply to look at the x^, 
defined as 



X 



(24) 
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TABLE L The amplitudes and their significance for different methods of calculating the covariance. The left columns show 
the two Monte Carlo methods, while the right two show the jack-knife method with equal area (10 deg^), in one case masking 
only patches of the density map, and in the other masking both density and temperature maps. We do not calculate the full 
covariance matrix or the total significance for the jack-knife cases. For 2MASS, we have cut the first four angular bins because 
of their SZ contamination; the total significance is obtained discarding these bins. 



where the inverse covariance matrix and the data can be 
referred either to a single catalogue or to the total mea- 
sure. Whereas the likelihood method discussed above 
looks at how well a model can reduce ? it is also worth 
simply looking at the magnitude of for the null hy- 
pothesis test, where we calculate the Xo assuming the 
theoretical cross-correlation is zero. 



catalogue 


/ 


xi 


•1 

Xbcstfit 


'2 

Xacdm 


2MASS 


9 


5.4 


5.2 


5.2 


SDSS 


13 


17 


11 


12 


LRG 


13 


9.6 


4.9 


5.7 


NVSS 


13 


17 


6.0 


6.3 


HEAO 


13 


18 


10 


10 


QSO 


13 


9.7 


3.7 


4.0 


TOTAL 


74 


67 


47 


48 



TABLE II; A comparison of the absolute for the various 
experiments. 



In Table [TTl we the show the x^ for the null hypothesis, 
as well as for the ACDM and best-fit models. We use the 
MC2 errors, dropping the first four bins of 2MASS which 
appear to be affected by SZ. While there is much varia- 
tion, in most cases there is not clear evidence against the 
null hypothesis, in that its Xo is not significantly greater 
than the number of data points. However, the x^ values 
are significantly reduced if one assumes one of the models, 
like ACDM , which predict a non-zero cross-correlation. 

The reasons for the particularly low x^ for the LRG 
case is unclear, and we investigate this more below. It 
might be an indication that the error estimates are in 
some sense too large, or that the covariance between 
angular bins is different than expected from the sim- 
ple Monte Carlo simulations, perhaps as a result of fore- 
grounds. However, it should be emphasised that the x^ 
for the null hypothesis is fairly conservative, and unlike 
the Bayesian likelihood approach, it fails to account for 
the fact that we have strong theoretical expectations for 
the signal we are looking for. 



D. Eigenmode decomposition 

To better understand the covariance of our data, and 
especially to understand the x^, it is useful to study the 
eigenmode decomposition of the covariance matrix. As 
a worst-case example, we will use here the measurement 
and covariance matrix for the LRG sample calculated 
with the MC2 method (dimension n = 13). 

We can factorise the covariance matrix into the form 

n 

C^, = E ^l^^i^^^ (25) 

where Ay = bij is a diagonal matrix whose elements are 
the eigenvalues of ; the rows of J7y are the 13 eigenvec- 
tors of the covariance matrix. We plot the variances, 
Ai = cr|, in the top panel of Fig. [71 and some of the eigen- 
vectors are shown in the bottom panel. There, we can see 
that the modes associated with the biggest variance are 
the low frequency ones, while the low variance modes os- 
cillate significantly. This reflects the fact that the great- 
est differences between the Monte Carlo realisations is 
in the low frequency behaviour of the cross-correlation 
functions. 

Both the measured and theoretical CCFs can be de- 
composed into this eigenvector basis. In particular, any 
cross-correlation vector can be written as v = Ajej, 
where Ai = v • ej. We show in Fig. [5] the decomposi- 
tion of the data and theory divided by the square root 
of the variance, <Ji. For a typical CCF from the A/fonte 
Carlos, these amplitudes should be Gaussian distributed 
with unit variance. We can see how the smooth shape of 
the theoretical real space CCF is reflected in this eigen- 
mode decomposition: the theoretical amplitude is very 
well approximated by the first two modes only. However, 
this is not the case for the measured CCF, for which 
higher frequency modes are also significant. 

We next look at the contributions to the x^ from the 
different eigenmodes. We show in Fig. [9] the evolution of 
the cumulative xf , i-e. the cumulative contribution to the 
X^ from each eigenmode. Here we compare the raw x^ 
from the observed cross-correlation function to that for 
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FIG. 7: Eigenvalues of the MC2 covariance matrices of the 
cross- correlation between the LRG sample and the CMB (top 
panel), and first three eigenvectors (bottom panel). The red 
dashed line shows the highest frequency mode. 



the residuals when the theoretical models (ACDM and 
the best fit amplitude) are subtracted off. As expected, 
the theoretical models only impact the lowest two eigen- 
modes. The low x^, however, is largely the result of 
the higher frequency modes, which seem to have slightly 
lower amplitudes than is seen in the Monte Carlos. 

If we consider only those two modes which are expected 
theoretically, the for the null hypothesis is actually 



fairly high: X2 = 4.8. This would exclude the null hy- 
pothesis at more than the 90% level. 
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FIG. 8: Eigenmode decomposition of the amplitude of the 
measured (red dashed), theoretical (black solid) and best fit 
(green long dashed) CCF. 



VI. COSMOLOGICAL CONSTRAINTS 

Assuming the observed cross-correlations are produced 
by the ISW effect, we can compare them with the the- 
ory predictions to obtain cosmological constraints. As 
described above, the ISW temperature anisotropics are 
produced as a result of time variation in the gravitational 
potential, and it is the evolution of the potential which 
our measurements constrain most directly. The cosmo- 
logical parameters which impact the linear evolution of 
the potential are the dark energy density and its evolu- 
tion, and the curvature of the Universe. 

The actual cross-correlation measurements will also de- 
pend on the nature of the large-scale structure probe, its 
spectrum and its bias. For example, if we normalise to 
the large scale CMB, changing the shape of the power 
spectrum (e.g., by changing the Hubble constant or the 
dark matter density) will change the variance of the dark 
matter distribution on smaller scales, quantified by erg. 
Since the ACFs of the surveys are fixed by observations, 
changing trg effectively means a different bias will be in- 
ferred for each survey. 

The cross-correlations will rise and fall with the 
amount of structure in the probe. Thus, instead one 
could focus on the dimensionless cross-correlation, effec- 
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FIG. 9: Cumulative Xi obtained summing the contribution up 
to the i-th eigenmode, for the three models: null hypothesis 
(red solid), best fit and ACDM . 



late-time ISW effect. Our observations of the ISW cross- 
correlations rule out such models at the ~ 4.5(t level, 
based on the difference in the between the null hy- 
pothesis and the ACDM model in Table HIl Such models 
also struggle to fit the recent observations of the angular 
scale of baryon oscillations (b^] . 



B. Flat ACDM models 

Next, we study the likelihood of a family of flat mod- 
els with varying 17™, f^A = 1 ^ ^m- As we can see in 
Fig. [TUl the ACDM model is an excellent fit to our data: 
the IfT interval for the parameter is Sim = O-^S^oor us- 
ing the MCI covariance estimate. A higher ISW signal 
and slightly lower estimate for Sl^ results from the MC2 
errors (Sim — 0.20lg;Q7); this is due to the higher best 
fit amplitude in this case. The error bars can be seen 
to be very asymmetric, as the ISW effect increases dra- 
matically when the matter density becomes small. Mod- 
els with S7„j < 0.1 would predict a much greater cross- 
correlation than is observed. 



tively 



Ir^TT/^gg I r'TTr'SS ' 

which removes the dependence on bias (assuming it is 
linear) and probes more directly the ISW effect itself. 
The ISW effect arises on fairly large scales, e.g. k ~ 
0.01ft.Mpc~^, depending on the redshift distribution of 
the survey. Equivalently, for each model, we calculate 
the bias of each survey based on the observed ACF, and 
use this to find the predicted CCF for the model. 

This makes our measurements largely independent of 
parameters other than fli:iE,w,Cs and flk- In practise, 
we choose to keep the dark matter physical density fixed 
uj,n = Slm/i^ = 0.128 to the WMAP best fit value, but the 
constraints are largely independent of this assumption. 



A. Models without dark energy 

While many independent probes seem to indicate the 
existence of dark energy, it is worth exploring models 
which might account for the observations without dark 
energy; recently, an attempt has been made that does 
this, but which requires a significantly lower Hubble 
constant, modifications to the primordial power spec- 
trum and other non-standard features |6l|. Such mod- 
els would be dark matter dominated today, and have no 



C. Flat wCDM models 

We next study the likelihood of a family of flat dark 
energy models, where we allow the dark energy density 
to evolve with equation of state w. The results are shown 
in Fig. [Til from which we can see that ACDM {w — —1) 
is very consistent with the measures. We can understand 
this if we observe that the measured excess in the ISW 
signal is largely redshift independent, while models along 
the same degeneracy line with a lower (higher) w would 
predict an excess at low (high) redshifts respectively. 

Initially we assume the dark energy sound speed is = 
1, as is typical in scalar field models like quintessence. We 
also show the same range of models, but with a different 
dark energy sound speed = in Fig. [11] We can see 
that in this case the degeneracy line changes direction due 
to the clustering of dark energy. ACDM is still a good fit 
to the data, as the cosmological constant likelihoods are 
not affected by the sound speed, and there is no clustering 
in that case. 

The constraint on the sound speed itself is very weak. 
There are too many dark energy parameters (density, 
equation of state, sound speed) to expect any strong con- 
straint. We reduce the numbers by assuming the CMB 
shift parameter is fixed to the observed value, coupling 
the equation of state to the dark energy density. The 
results can be seen in Fig. [121 Even with this additional 
constraint, the sound speed is weakly constrained be- 
cause the data are consistent with a ACDM model, where 
there is no dependence on the sound speed possible. 
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FIG. 10: Likelihood for flat models with varying Qrn from the 
MCI and MC2 errors. The shaded areas represent 1, 2 and 3 
a intervals for Qrn- ACDM is a good fit to the data. 



D. Curved ACDM models 

Since curvature can also cause the gravitational po- 
tential to evolve, we explore the constraints if we relax 
the flatness condition. However, for simplicity we as- 
sume the dark energy is a cosmological constant. We 
study the likelihood of fim, with a corresponding curva- 
ture f2fc = 1 — r^A — ^rn- We explore the full Qm ~ ^^A 
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FIG. 11: Likelihood for flat models with varying flm and w 
from the MC2 errors. The shaded areas represent 1 and 2 
a intervals. The top panel assumes relativistic sound speed, 
such as would occur in a quintessence model, while the lower 
panel assumes the opposite extreme of zero sound speed. 



space; we see the relative likelihoods in Fig. [T3l which is 
obtained with MC2 errors. 

From this figure, we see that ACDM is still a good fit 
to the data. An interesting feature of this figure is the 
degeneracy line between Vim and this is related to 
the relative efficiency of the curvature and dark energy 
as sources of ISW. Closed models (above the flat line) 
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FIG. 12: Likelihood for flat models with dynamical dark en- 
ergy as a function of the sound speed, where we fix the matter 
density based on the equation of state, assuming the CMB 
shift constraint. 1 and 2 a intervals are shown. No constraint 
is possible for the cosmological constant limit {w = —1). 



FIG. 13: Likelihood for curved models with varying Qm and 
57a from the MC2 errors. The shaded areas represent 1 and 
2 a intervals. ACDM is a good fit to the data. 



give negative ISW, and can cancel the effect of increasing 
the cosmological constant, while the opposite happens for 
open models (below the flat line). 

E. Comparison with other constraints 

Finally, we wish to compare the ISW constraints to 
those arising from other cosmological observations, in- 
cluding the CMB power spectrum, baryon oscillations 
and type la supernovae. For the latter, we use mea- 
surements of the luminosity distance from the Supernova 
Legacy Survey [Ij. 

For the CMB observations, most of the dark energy 
information (at least that independent of the ISW effect) 
is distilled in the CMB shift parameter, defined as 

R=^nZHo-{l + z,)dA{z.), (27) 

where dA{z) is the angular diameter distance and is 
the redshift of the last scattering surface {z^, = 1090); 
this expression in the flat case reduces to 

R^Vn7.Hol (28) 

R has been measured to be i? = 1.70 ± 0.03 We can 
see from Fig. [T3] that this constraint has a degeneracy 
direction parallel to the ISW degeneracy in the flat case, 
but is less so in the general curved case. 



Finally, for the baryon oscillation (BAO) measure- 
ments [0] we use the constraint on the volume distance 
measure defined as 

dviz) = [(1 + z)24(z)z c/i/(z)]i/3, (29) 

The constraint on this parameter by Q is 
dy(0.35)/dy(0.2) = 1.812± 0.060. 

The SN data is orthogonal to the ISW constraints, 
and jointly they are consistent with the ACDM model; 
there is little evidence for additional curvature or evolv- 
ing dark energy. The CMB shift constraint is similarly 
consistent with the cosmological constant concordance 
model, though the constraints are not as orthogonal to 
the ISW constraints. The ACDM model preferred by the 
SN and ISW measurements is consistent with the CMB 
shift combined with the measurements of the Hubble con- 
stant from the HST Key Project (5^]. 

The exception to this concordance picture comes when 
the BAO data is considered. The BAO contours are 
similar to those from the SN, but shifted. In the flat 
dark energy case, the combination with the ISW prefers 
a larger dark energy density which has increased with 
time (phantom). When all observations are combined, 
the BAO data are swamped by the SN data, and the 
result is fully consistent with the concordance model as 
found by Q. 
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FIG. 14: Comparison with constraints from other observa- 
tions, including CMB shift (black), SNe (red) and BAO (blue) 
(top panel), and combined likelihoods using the ISW + each 
one of these other constraints (bottom panel, same colour cod- 
ing). 1 and 2 a contours are shown (solid and dashed lines 
respectively). The MC2 errors are used. 



VII. CONCLUSIONS 

In this paper, we have measured the cross-correlation 
between the CMB and a large range of probes of the 
density in a consistent way, and have calculated their 
covariance taking into account their overlapping sky cov- 
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FIG. 15: Same as Fig. [T4]for the curved case. 



erage and redshift distributions. While individual mea- 
surements vary somewhat depending on how the data are 
cleaned and how the covariance is calculated, the overall 
significance of the detection of cross-correlations is at the 
~ 4.5(7 level. 

These observations provide important independent ev- 
idence for the existence and nature of the dark energy. 
The observed cross-correlations are consistent with the 
expected signal arising from the integrated Sachs- Wolfe 
effect in the concordance model with a cosmological con- 



stant. The observed signal is slightly higher than ex- 
pected, higher than the expectation from WMAP best fit 
model by about Icr, thus favouring models with a lower 
f2m- However, we do not see any significant trend for the 
excess as a function of redshift, and so there is no indi- 
cation of an evolving dark energy density. By combining 
these results with other cosmological data, we find a gen- 
erally consistent picture of the behaviour of the Universe, 
which is converging towards the ACDM model although 
the uncertainties remain considerable. The only partial 
exception to this picture is the BAO result which, even 
when combined with our ISW measurement, is in slight 
tension with the ACDM model (at ~ la). 

The results of our analysis and the covariance matri- 
ces are available upon request from the authors (contact 
tommaso . giainnEintonioQport . ac.uk). 
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APPENDIX A: CORRELATED MONTE CARLO 
MAPS 

1. Basics 

Here we describe how to make Gaussian maps with 
a prescribed set of auto- and cross-correlation functions 
for use in the estimation of covariance matrices. Let us 
assume we have n maps, which could include tempera- 
ture and various density maps at different redshifts or 
frequencies. Let us call these maps nii where i ranges 
from 1 to n. 

Any two maps, nii and rcij , will be correlated and these 
correlations will be described by a correlation function 
C"^ {i}) and associated multipole moments . These 
correlations will be symmetric under interchange of the 
maps, C"-' (7?) = C^' (■(?), so we have n(n+l)/2 correlation 
functions or spectra which describe the two maps. 

Most map making algorithms, like synf ast[37j. work 
in Fourier or spherical harmonic space. Effectively every 
mode is given a random amplitude ^, which is a complex 
number with unit variance and zero mean: {£,£,*) — 1 
and (^) = 0. These are then multiplied by the square 
root of the power spectrum in order to ensure the proper 
correlation functions. (There are additional constraints 
to preserve the reality of the fields on the lattice, e.g. 

= ^^jj., but it is not necessary to go through these 
here) . 

It is sufficient to consider a single mode or harmonic 
amplitude of each map, as all the others will be similar 
but independent. Assuming we are working with spheri- 
cal harmonics, we want to ensure that 



(Al) 



The S functions follow simply from using uncorrelated 
random amplitudes for each harmonic mode. For a single 
map, the right power spectrum is ensured by simply using 

(A2) 

and this is effectively the prescription used by synfast. 

When considering more maps, it is necessary to use 
more random phases, building the final maps from a 
combination of different maps. With n maps, n diff'er- 
ent phases are required for each mode. Here, we de- 
note the different phases with Latin letters, a, 6,c, ... 
Different phases will be assumed uncorrelated, so that 

The simplest example is to consider two correlated 
maps, mi and ni2. These are described by three spectra: 
C|^, Cj^ and Cf^. These are made using the amplitudes 



It is simple to verify that with these amplitudes. 



This is simple to implement with synfast. First cre- 
ate a map with power spectrum C|^, and then make a 
second map using the same seeds and power spectrum 
(C|^)^/C|^. Add this second map to a third map made 
with a new seed and with power — {C^)'^ / C}^ . Note 
that this should never require taking the square root of a 
negative number; however, if its very strongly correlated, 
numerical errors could cause problems. However, for the 
weak correlations considered here, this is never an issue. 

The only difficulty is that this inherently produces pos- 
itive correlations, as the default of the synfast code. 
This can be worked around simply. For example, if 
is always negative, one can simply flip the signs of the 
second map after it is produced. If instead changes 
sign, then break up the power spectrum into positive and 
negative pieces, making a map for each and subtracting 
the 'negative' map from the 'positive' map. 



2. The general case 

Next we consider an arbitrary number of maps. For 
simplicity, we drop the £ and m subscripts where the 
meaning is unambiguous. Effectively, the challenge is to 
solve for a particular set of amplitudes T, where 

= £,aTla 

— ^aT2a + £,bT2b 

3 



a 



a 



^aTsa + £,bT3b + ^.cT^c 
^aTia + ^bTib + ^cTic + ^dTid (A4) 



etc., subject to the constraints that (a^a^*) = C"^. 

One thus has n{n -I- 1)/2 equations with the same num- 
ber of unknowns T. These begin as: 



^22 

^23 
^33 
^14 
^24 
^34 
^44 



_ rp2 

= TiaT2a 

r7-T2 

— -'2a 
= TiaT'g.a 

— T2aTza 

= Tl+T, 



Tib 



T2hT, 



3b 



TlaTia 
T2aTia 



26-t3h 
T2bT4b 

T^bTib 



T, 



4b 



Tf 



id 



(A5) 



etc. While quadratic, these can be solve d in stages lin- 
early. Solve the first for Ti g = Vc^. Use the sec- 
ond to show, T2a — C^^/VC^i and the third to get 
T2b = V'C22 _ (c'i2)2/C'". This reproduces what was 
shown above. 

After this, things continue similarly. At each point, 
we use the next equation to solve for the next missing 
variable: 

Tla = VC^ 



T2b = ^C'22-(C12)2/C11 

T36 = ((723 _ Cl^C^/C")/ VC22 _ (^12)2/^11 

//o23 /^12^13 



1/2 



^3c / / ^ ^22 _ (^(712-)2 

Tib = (C24 _ C^^C^^/C^^)/^C^^ - (Ci2)2/(7ii (A6) 

etc. Things will take similar forms as one goes on, but 
getting progressively more complicated. 

It can also be programmed recursively, which may be 
simpler to implement. By this, we mean. 









T2a 




C^VTla 


T2b 












Tsb 




(C^^ _ T2aT3a)/T2b 












C^^/T^a 


Tib 




((^24 _ T2aTia)/T2b 


Tic 




((734 _ T^^Tia — T3hTib)/T3c 


Tid 




Jc^^-Tl-Tl-Tl 



(A7) 



etc., with each step using only variables already solved. 
The general recursive expression for these spectra is 



T- ■ — 



T- ■ = 



k=l 

c^'-i:rJ,T,kTjk 



if i > j. (A8) 



"-33 



These amplitudes arc squared for the input spec- 
tra for synfast, but one must beware negative cross- 
correlations as discussed above. A simple modification 
to a program like synfast could enable it to read in 
amplitudes rather than spectra, and this would be more 
efficient compared to reversing the sign of the maps after 
they are created. 
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